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SUMMARY 

The point-souroe-diatrihution. method of calculating the 
aerodynamic ooeff Icler^ts of thin vings at supersonic speeds was 
extended to include the effect of the region hetwepn the wing ^ 
hound£iry and the foremost Mach waye from the wing leading edge. 

The eff.ect of this region on the surface velocity potential has 
been dete\TiilnBd by an equivalent function, which is evaluated 
over a portion of the wing surface. In this manner, the effect 
of angles of attack and yaw as well as the asymmetry of top and 
bottom wing surfaces may be cdloulated. As examples of the method 
the pressure distribution on a thin plate wing of rectangular plan 
fora as well as the lift and the drag coefficients as a function 
of -Mach mmber, angle of attack, and aspect ratio are calculated. 
The equations for the surface velocity potential of several other 
plan forms are also included. 


INTRODUCTION 

The theoi’etical and experimental investigations of aircraft 
performance at supersonic speed.3 have been greatly stimulated by 
modem developments in high-speed flight. The theoretical aero- 
dynamic perfonnanoe of thin wings nevertheless has not been com- 
pletely solved, even throu^ the approximations of the line arized 
Prandtl-Glauert equation. 

Puckett (reference l), by means of a point source distribu- 
tion, has formulated a method to derive the pressure distribution, 
the wave lift, and the w'ave drag for thin wings at angle of , 
attack, pmvlded that the leading edge or the wing tip, as the 
case may be, is swept ahead of the Mach line . The method gener- 
ally falls when the sweepbaok is greater than the Mach line 
because the flow over one surface of the wing can Isafluence the 
flow on the other surface. 
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Jones (reference 2) has been able to calculate the pressure 
distributions on a series of supersonic wings by means of line 
sources. The results are the same as would be obtained by 
Puckett's theory, however, and the method is subject to the same 
limitations . 

After transforming the Prandtl-Glauert equation to curvi- 
linear coordinates, Stewart (reference 3) picked the special 
solutions corresponding to conioal flows. In this manner, the 
Prandtl-Glauert equation was reduced to the~T;wo -dimensional 
Laplace equation that permits the use of conformal mapping. As 
a special case, Stewart obtained the lift distribution on a thin 
delta wing at small angles of attack. Brown (reference 4) has 
independently solved the same problem by use of a doublet line 
source distribution on the wing surface. 

The present paper extends the point -source -distribution 
method (applied by Puckett to the wing surface) to include the 
effect of the region between the wing boundary and the foremost 
Mach wave from the leading edge. By use of a source distribution 
external to the wing, the intoi’action of the two wing surfaces 
may be isolated. In this manner the pressure distribution in the 
vicinity of the wing tip, as well as the effect of profile shape, 
angles of attack and yaw, and aspect ratio, may be calculated for 
a series of finite wings. This work was performed during Febniary 
1947 at the NACA Cleveland laboratory. 


ANALYSIS OF MEIEOD 

Thin wings will be so used in the analysis that, the perturba- 
tion velocity components may be assumed to be small compared to 
the free-stream velocity. The linearized partial differential 
equation for the velocity potential of a ccmprossible fluid may 
then be applied. The problem is to find a perturbation velocity 
potential that will; (a) satisfy the linearized partial differ- 
ential equation of the flow, (b) vanish in the region ahead of 
the foremost Mach wave, (c) give streamlines that are tangent to 
the airfoil surfaces, and (d) take into account the interaction 
between the top and bottom wing surfaces as represented by the 
perturbed field between the wing boundary and the foremost Mach 
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The Prarn‘’.tl- G-lau vt linearized equation for the velocity 
potential of a nonvl 3 < o /‘3 irrctational c cmprese lb le fluid may be 
vritten as 

,1 ^ = 0 ( 1 ) 


where 

M f roe-stream M acVi mmDer (undistiurbed flow parallel 

to x-axia) 

cp pGrt^U’ba^ion velocity potential 

X or I'l 

y or T] > Cartesian coordinates 
z or ^ I 

For convenience the sjyribols ar-; defined in appendix A. A basic 
solution for the potential, of a unit point source disturbance at 
U is . 


~l (2) 

^ /'U - 1)2 - ^2 (y - 32 (z - t )2 


•where 

3 = “ 1 

More gr-.j.eral eolut.ions my be obtained by integi’ation to givo 




g»d| dT]dt 


Jj 


(3) 


w*'.ere q.' is th,c source strength por unit volume. For the -thin 
vlng, the soui’ce>£ and the wing may lie in the x, y plane and thus 
4 'd^ may be rejpl^ced by q., the secures strength per xmit area. 
Eq.uation (3.) then becomes 
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Cp 




A 


X 


~X)‘ 


■M-ln ^ - 

- 3“ (y - 


( 4 ) 


Puckebt (reference l) has shown that the houndary conditions for 
thin wings may he. satisfied as z -» 0 by setting q. = w/rt, where 
w is the perturbation velocity component normal to the x, y plane. 
The quantity w is proportional- to the loc^l slope (the angle sub- 
tended by the wing'' surface frcaa the x, y plane in tj = constant 
planes) of the wing in the free-atream direction at the point (i>Ti). 
If 'h represents this slope, - 



(5) 


where U is the free-stream velocity. Equation (4) then becomes 


cp 


U 

It 


Xdt dn 


y (X - if - 3My - nV 


( 6 ) 


The form and the derivation of equation (6) indicates that 
alteration of the local slope at point (I 2.1 will not 

change the perturbation velocity component w at some other 
point (|, p). The velocity potential at any point (x, y) on 
the surface of the wing may then be calculated by integrating 
equation (6) over the region in the x, y plane bovinded by the 
forward Mach cones. (Bee, for example, fig. 1(b).) Puckett 
restricted hie integration to the wing surface, where X is 
assumed known. The solutions obtained in this manner are valid 
if the wing is swept less than the Mach angle or if the top and 
bottom surfaces of the wing for any sweepbaok angle eire symmetrical 
about the x, y plane. - 


If a proper distribution of source strength XU/it is chosen 
for the regions between the foromost Mach wave and the leading 
edge, equation (6) will give the velocity potential at any point 
(x, y, 0) regardless of sweepback angle and asymmetry of top and 
bottom wing surfaces. 


The strength of the source distribution between the Mach 
cone and the leading edge (or wing tip) must correspond to the 



NACA TN Wo . 1382 


5 


local perturbat j on velocity component w of the region. This 
velocity is in turn influenced .ly the slope of both the top and 
the bottom surfaces of the ving. 

A thin impermeable diaphragm is assumed to coincide with 
a stream sheet in the x, j plane between the wing surface and 
the foremost Mach wave. The presence of the diaphragm will not 
alter the flow over the wing surface i The diaphragm may then 
be regarded as an extension of the wing to eliminate the 
external field between the wing boundery and the foremost Mach 
wave. 


Because the diaphragm coincides with a stream sheet, it may 
sustain no pressure difference at any point between its top and 
boutom surfaces. Furthermore, there can be no discontinuity in 
the velocity components across the diaphraga. This situation 
• requires that the surface velocity potential at any point on the 
top and bottom surfaces of the diajhragn. are eq^ual. Inasmuch 
as the extended wing allows no interaction between its two sur- 
faces, the velocity potential at any point (x, y) may be 
calculated from eithei- the top or the bottcan auif ace of the 
original wing and diaphiagm. 

The local slopes of the wing on its top and bottom sur- 
faces at the point ( |, t]) may be represented by and Og, 

and A may represent the corresponding slope on the top sur- 
face of the diaphragm. (For’ oonyenience, the sign of a is 
oppositely defined on the two surfaces. For example, aip and 

Og are both positive on a wedge profile at an angle of attack 

of 0.) The areeis of the wing emd the diaphragm surfaces 
included in the forwai’d Mach cone frcn a point on either the 
wing or the diaphragm are represented as S-^ and Sp, respec- 
tively. Number subscripts 1, 2, • * • represent sections of 
each of these areas. Then hy equation (6), the velocity 
potential at (xp^yp) (fig. 1(a) ) is 



cp 




U 

ir 


'i’l 


OJ 


Sy V^(XD - if - P“(JD - n)“ 


u 

rt 


■^(|y_n) ^_4 il 


Sd v(xi> - i f - p“(yj3 - n)' 




U 


aB^-E 


n // »■ \2 „ 2 /_ 

°v V V-»-D - Zf - P Ub “ ’1^ 


_ _ U pp 

" L /C 


~\(t, Ti) at '^n 


c\2 „2/^ ^2 

JJ^ - P 


viiei'e T>d^t ^,B potentials on the diaphragm calculated from the top and the 

■bottom surfaces of the wing and the diaphragn, respectively. Or 


‘ ?^L.JliJj.-.^Ii 

Sp -if- 3^(71, - J 


( ob - <7t) 

Sy 2 y(iB - if - (7d 7 


(7) 


This int.egrfli e'nuation defines the function X, The velocity potential at any point (x^ y) 
on the top surface of the wing is given from equation (6) as 


9.--U 

J, 


Oipd£ dti 


— It 

tSf V “ P"^y - Tr 


u '■'<■' 

ir 


Hi, Ti) 


. _ ! t *\2 

j aj) V tx - ir - \i~Kj ~ T)j“ 


( 8 ) 


similarly, the potential on tdie bottom surface is 




u 

rt 


rr. ^ t Am 


+ li 




7(x - if - p2(y _ 1^7^ j 


^ y * % ^ #r -I 


Sd 


y(i - ^ 


8^(7 - 


(8a) 
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vhere the integration is made over the region in the y plane hounded by the forwaini Mach 
cone frail the point (r, y) over ail wing and diaphragm surfaces. 

In special cases, the potential op of equation (8) may b e obtained wit hout explicitly 
solving for This simplification is- most easily accomplished iu an oblique u^ v 

coordinate system the axes of Tdilch lie parallel to the Mach waves. In this systotj, for , 
example, the value a£ the coordinate u ia the distance frcm the v-axls to the point 
measured parallel to the u-axis. The transformation equations are 


e - |(t + u) 

Inasmuch as the elemental area in the u, v 
be written either as 


(e + 3Tl) 

11 = ^ (r - u) (9) 

20 

coordinate system is du dv, equation (7) may 


20 


Xfu, v) du dv 20 


(cr-R - aj) du dv 

2 

M 


sj, ^ 


3^ 


(7a) 


'j 


or as 
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M 




y) du dY 

J{^D - n) (tj) - 


1 

M 




(C^B - a^i) du dv 
2 aA'^D - Ti) (Vj) - v) 


(7b) 


where Uj^ and represent the coordinates of the point (Xp, y^) In the vl , y systeua. 

, j 

The regions of integration and the coordinate systems of equation (7) are ajcetched for a 
wing plan riew In figure l(a). The zero of the coordinate system is placed at the point of 
tangency of the foremost Mach wave and the leading edge. The wing area is bounded by the 

two curves (or two branches of the aaine curve) v = v^ (u) and v = Vp (u) and the 

M ^ ^ 

line u = Up = ( 2 p - 3yp). Application of equation 7(b) to this case yields 


CD 


aa- 

?v(u, v) dv ^ 

ru-n 

® du 

1 (Gp -Gqi) av 

- li) 

Jo ^ 

v-2(u) 

«/(Up - u) 

w 

CO 

! 

r 1 
' ^ 


( 10 ) 


Inasmuch as the limits of Integration of the 


du 


integrals are same for all values ■ 


of Up and orwitig to the nature of the functions, the two integrations with respect to v may 
be ©quat^ along lines of constant Vp that extend across the wing and the diaphragn. 




V2(u) 


X(u, v) dv 


r^ 2 (u) 

■^l{u) 


(gp - Cip) dv 


( 11 ) 
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!33ie contribution to the Telocity potential on the top 'surface of 
the wing attributed to the diaphragm, (fig. 1(b)) is given by 
eq^uations (8) and (9) as 




U_ 

Mjt 


■u' 


du 


y(u^ - u) 


'w 


Tp(u) 


-/(Tv - t) 


( 12 ) 


^ere u^ and 
wing and the limit 


are the coordinates of point (x, y) on the 
u' is obtained by solving the equation 


V^ = Tg (U* ) 


Ihe integration limits with respect to v and the Integrand of 
equation (12) are the same as the left side cf equation (11), 
except that v.^ replaces Vj^; but the value of Vjj along the 

T = constant line peiseing througjh the point (u^, v^) is v^.. 

The second member of equation (11) may therefore be substituted 
into equation (12) to give 


'^^T,D = “ 


Mjt 


nu.' 


du 


0 - u) 


Pn - °T) dv 

•®'l(u) 


^ yc'^w - v) 


It 


(C^ -Oiji) dg dq 


1^,2 2 




tf - 3^(y - nf 


( 13 ) 


The contribution of the diaphragm to the potential on the 
wing surface may thus be replaced by an equivalent integration 
over a portion of the wing surfaoe. Ihe potential on the wing 
surfaoe is then 
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cPgn = 


U 

ir 




^v ,2 2 a/^ - i)^ - 3^(7 - r\)^ 


or 


U 

jt 


■P 


Jm dE dTl 


K 

J. 


Sw,2 - E)^ - P^(y - I'l)^ 


Orp i^E 

V,;, -/(X - E)^ - P^(y - n)‘ 


u 

It 


u 

It 


O fp dE dll 


^ix - if - p^(y - ti)^ 


(og H-ax) (14) 


2 a/ (x - E)^ - p^(y - Ti)^ 


= 


2311 

itM^ 




du dv 




' (Or -f gy) du dv ^4^) 

Sv,2 - |(v + u)]2 " [y - 

The derivation of equation (13) includes the assumptionn of the 
linearized theory and the assumption that the leading edge is not 
blunt (corresponding to the use of a thin diaphragm). Aside from 
those restrictions j the equation includes the effect of asymmetry 
between the top and the bottom wing surfaces;. It may therefore bo 
applied to determine wave lift, drag, and pressure distribution in 
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the vicinity of wing tips of fairly general chordwise slope dis- 
trihutions. Because the only restriction on the functions (u) 

and T 2 (u) was that Sj) he influenced only hy the wing sec- 
tion Sy, the aerodynamic properties of fairly general plan forms 
may he evaluated. (In cases of so-called subsonic trailing edges, 
the solution for the velocity potential that is obtained vio- 
lates the Kutta-Jou!!row8kt condition in the vicinity of the trail- 
ing edge. The solutions may not correspond to actual flows under 
these conditions.) The affect of yawing the wing may also he 
determined simply hy simultaneously adjusting the functions V]_ (u), 

'’’2 *^T ^ amount corresponding to the angle of 

yaw. The effectiveness of wing tips and hence the effect of aspect 
ratio may likewise he detemined. 


V 




7 




\ 


t 


It 


EXAMPLES OF MEEEEOD 

Thin flat plate wing with rectangular plan form and no sweep 
hack. - For the flat plate wing (fig. 2), a-Q = - = angle of 

attack a and etiuation (14) heccmes 


'Pt =•— 


dn 

Swa-^Cs: - |)^ - (7 “ i))‘ 


(15) 


Thus, the external field Sjj cancels the effect of the region Sy^2 
as far as the potential at point (x, y) is concerned. 


The pressure coefficient in the region of the wing tip 

may he computed from the eguaticaa 



The value of Cp obtained from equations (15) and (16) is derived 
in appendix B for th^top surface of the wing to give 


m = - 3; + ^ sin 


-1 


^p,T 


0 * 



(17) 


t 
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(Eqtuation (17) is eq.uation (Bl) in appendix B.) The coordinate ,y 
iSj of ooTirse, negative. The preestire is therefore constant along 
radial lines from the origin^ has the free -stream value along the 
tipj and has the Aokeret value (reference 5) along the Mach line 
lying on the wing from the tip and leading, edge intersection. If 
the influence of the external field Sjj has "been neglected, the 

pressure coefficient would he one-half the Ackeret value along the 
wing tip instead of the correct value of 0. The result presented 
in eq.uation (17) was first derived hy Busemann (reference 6) and 
has been cited in reference 7. 


I 


The pressures on the top and bottom surfaces of the region 
influenced by the wing tip are integi'ated in appendix B to give 
the lift and dragboofficlents. The lift and drag coefficients 
are one -half the values obtained by the Ackeret theory (refer- 
ence 5). The wave lift and drag coefficients for the whole wing 
are given in terms of the aspect ratio A (if Ag-|) as 




Ct, = 


8 2PA J 


■'ll 


4a 

P 


(" - m) 


(18) 


which is derived as eq^uation (B4) in appendix B. This effect of 
aspect ratio on the thin flat plate wing has been previously 
reported in reference 8. 


Dlscontinuously swept wing of small finite thlokness except 
on edges . - The leading edge may lie on lines v = - k^u, and 

V = kgu, wiiere k-j_ and kg are positive constants. (See fig. 3.) 
For this case, equation (14a) becomes 


cPt = - 


2UP 


M 


du 


'kgu 


(g-R + gj) dv 


u)]2 - 1 (v-u)]' 


- 


'M 


du 


M 

‘28k. 


-(x+py) 


2P 


(x+py) 


g^ dy 




2 

( 19 ) 
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Similarly, the potential on the bottom of the wing is 


(g-R + ot) <^'7 


gB 

y^i (v+u)]® -p2jv -i (v-u)]2 

(19a) 

Equations (19) and (19a) apply for wings at angles of attack 
even thou^ the top and bottom surfaces are asymmetrical. For 
symmetrical profiles at an angle of attack of 0, they reduce to 
the expressions obtained by Puckett's theory. 

An interesting observation is that only the second of the two 
integrals in equations (19) and (19a) includes the effect of angle 
of attack j at an angle of attack a, ag = gg’ + a, g^ = aij' - a, 

and gg + g^ = gg' +g 5 i', -where gg' and gq>’ are the , local 

•wing slopes on the bottom and top surfaces at an angle of attack 
of 0. The first integrals of equations (19) and (19a) are iden- 
tical, Therefore, only -the second integral contributes to 

Both integrals must be considered when pressure distribution 
or drag coefficients are desired. 

As an example of the i^e of equation (19), the velocity 
potential for a wedge wing (fig. 4) of constant wedge angle 2g 
parallel to the free-stream direction has been calculated in 
appendix C as equation (G2). This potential for the top surface 
of the wing is 




' M 
3(3 kg 


(2C4-Py) 


'kgu 


du 


J-kgu 


2Up 




du 


M 

23kp 


K2:+Py) 




-kj^u 
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Ua J [(1 + ki) X + (1 - ki)j3y ‘ ^ /k^ (x - py) 

" V (x + py) 


ua } 


1 


[(kg - 1) X - (kg + 1) py] y(kg - 1) X - (k 2 + 1) py ^ 
a/ ^ 2 ykg (x - py) + v^x +' py | 


np 


j / (^1 + ^2) + Py) [(^2 - 1)^^ - (k 2 + 1) p y] 


(1 + ki) X + (1 - ki) py 

aAi 


tem 


-1 


jki [(k2 - 1) X - (k2 + 1) py] 
(k]^ i- k 2 ) (x + py) ^ 


( 20 )' 


The Telocity potential for the hottom surface may be obtained by 
replacing a by -a. Only the second brace of equation (20) influ- 
ences the lift of the wing. Pressure coefficients may be obtained 
by substituting equation (20) into equation (16). 

Wing influenced by two independent perturbed flow fields 
external to wing surface . - An external flow field is considered 
to be independent if it does not include an external flow field of 
unknown strength in its forward Mach cone. (See fig. 5.) If 


K (x, y, i, Ti) £ 


U 


- |)2 - p2 (y _ 

the velocity potential at point (x, y) will be 

rr 


cpr 


aqiKd| dq - 


^w 


^■^^Kd^ dq 


^D,l 


AgKdl dri 


J^D,2 


( 21 ) 
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By application of equation (15) 



"V 

dq = 


s 

^D, 1 



(Ob -aqi) K 


d| dTi 


w,l 

(ag - Ot) K 
2 

=w,2 


d| dq 


and 


(22) 



'X 

XgEdt dq = 

ft ft 

*• >. 

Sd,2 



(0jB - ap) E 


d| dq 




(0B -Ot) g ,, , 
2 


S. 


V,2 


Also 


(23) 



ft 

ft 

ft 

a^d| dq + 

^ r* 

V* V. 



> 

s ^ 

V, 1 



a^d| dq 


^j2 


-H 

ft 

ft 

a^d| dq + 

•» ft 

u 

V 

Sv,3 

■k. 


OipEdl dq 


( 24 ) 


Sutetitution of equations (22), (23), and (24) into (21) yields 
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CpT = “ 

•V 

(0B +Ot) K . 

2 d| dq - 

"V p 



®W, 1 V 



Sv,Z 


a^Kd-l dTj 



■' (<7-r + Om) K 

p 


— d|dTl - 



Sw,3 



ajEal dll (21a) 




Similar extensions may "be made for other Independently perturbed 
flow regions in the forward Mach cone. (The boundary of the 
shaded region in figure 5 gives the limit of validity of equa- 
tion (21a). The shaded region is influenced by external fields 
that are no longer independent . ) 

Because equations (7) and (8) are linear with respect to the 
local wing slope o, equations (14), (19), and (21) may be 
divided into two sets of integrals; the first set will depend on 
the angle of attack a, but will be independent of the wing 
slope o' at an angle of attack of 0; the other set will be 
independent of angle of attack but will depend on the slope o ' . 
For symmetrical profiles about the x, y plane (at an angle of 
attack of 0), the second set consists of the velocity potential 
for the given plan foim at an angle of attack of 0. The first 
se't represents the velocity potential of the thin flat plate wing. 
For symmetrical profiles at angle of attack, the aerodynamic 
coefficients for the wing may therefore be obtained by super- 
posing the solution at an angle of attack of 0 (calculated by 
the methods of Puckett (reference l) or Jones (reference 2)) and 
the solution at angle of attack of a thin flat plate wing of the 
same plan fom. From the form of equations (7) and (8), this 
superposition is apparently general. In this manner, for example, 
solutions for the symmetrical delta wing and the related airfoils 
of somewhat arbitrary chordwise thickness distribution may be 
obtained from the results of Stewart (reference 3) or Brown 
(reference 4) . 


DISCUSSION CfF METHOD 

The general equations (7) and (8) will satisfy the boundary 
conditions of thin wings at supersonic speeds. The lees general 
solution of equation (13) may be applied to calculate the 
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oontrilDution to tlie velocity potential on tlie wing of a single 
perturbed field between the ^rtng boundary and the foremost per- 
turbed Mnxjh cone; as illustrated by equation (21a)^ the method 
may be extended to include the effects of a multiplicity of 
independently perturbed external flow fields . The basic 
equation (8) from which the succeeding equations were derived 
incl.udee the effects of asymmetry about the x, y plane between 
the top and bottcm surfaces, such as would occur at angles of 
attack. The method may also be used to calculate the effects 
of yawing the wing. 


Fli^t Propulsion Research Laboratory, 

National Advisory Conjmi.ttee for Aeronautics, 
Cleveland, Ohio, May 27, 1947. ' 
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APPENDIX A 
SYMBOLS 

The following syrnhols are used in this report: 


A 


aspect ratio 

- 



drag coefficient 

• 

Cl 


lift coefficient 


"p 


pressure coefficient 


1^1^ 

^2 

. constants greater than zero 


M 


free-atream Mach number 


4 


disturbance source strength per unit area 

* 

<!' 


disturbance source strength per unit volume 


S 


plan fom area 

m 

U 


free-stream velocity 


= 

V = 


\ obliq.ue coordinates whose axes lie parallel to 
Mach lines 

.. 

w 


z component of perturbation velocity 

- 

X, 


Cartesian coordinates 

. 



wing chord 


a 


angle of attack 

... 

P 


cotangent of free-stream Mach angle, \/m^ - 1 


t f 

t 

Cartesian coordinates 
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X 

a 


o' 

T 

Sutsoripts ; 

B 

T 

D 

w 


1, 2, 3 
Examples : 

Oiji 


Sw,3 

^1 

\ 

°P,T 


s lope of stream sheet near ^ ^ r\ plane 
measured in t] = constant planes 

slope of the wing surface with respect to the 
T) plane measured in tj = constant planes 

slope of wing surface at zero angle of attack 

perturbation velocity potaatial 

bottom (wing or diajhragm surface) 
top (wing or diajhragm surface) 
diaphragm (with exception of Cjj) 
wing 

refer either to numbered areas or curves 

slope on top wing surface 

potential on top surface of wing due to 
diaphragm 

wing area 3 

curve V = v^ (u) 

slope of diaphragm in plan area 1 

pressure coefficient on top surface of wing 
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APPENDIX B 

LIFT DISTRIBUTION NEAR TIP OF THIN WING 
OF RECTANGUTJ® PLAN FORM 

The llmita of Integrat/ion of equation (15) are evident from 
the following sketch. 



The potential at point (x, y) is then 
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Partial differentiation with r^pect to x yields 
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Therefore from equation (16) 

C = .?^=-5 + ^ sin"^ + l) 

P ■ rr riTF fl JtB ■ N. ^ . y 


(Bl) 
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The average lift coefficient 
for the flat plate wing as 


for .the wing tip may he 


computed 


Ct = 


^ 


2 dS 


Jas 


Because Cp is constant along rgidial lines from the origin, a 
triangular infinitesimal area is convenient. In terms of 


o 



the chord x. this area is 


dS 


Xq dy 


2 



4 



On the average, the wing tip area is one -half as effective aa the rest of the wing, although 
the lift-drag ratio in the frictionleas case is the same. In terms of the aspect ratio A, 
the wave lift and drag coefficients of the \diole wing are 


fi-X\ 

p V 2pAy 



(B4) 
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APPEND.TX C 

CALCULATION OF VELOCITI POTINTIAI. OF A DLSCONTINUOTJSLY 

SWEPT WEDGE WING 

For the wing of wedge angle 2a shown in figure 4_, equa- 
tion (19) hecoBies 
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Each of these integrals may "be integrated (reference 9, 
integrals 111 and 113), althou^ care should he teiten in the 
choice of signs for the stiuare roots. The following scheme 
was applied : 



\-h \ |-h' I = Iv'h >1'! 

' I - I I - I » ' - - I 

-v/ l\ = - Ia/i "h'l 

■ X i”t/\ »/ -I 

'\/-(-|b|)|‘b'‘l = 


iflaere d >^nfi d* are arbitrary mimbera. In otiier vorde, wien tvo negative slgie are multi- 
plied under the radical, the negative sign, ia transferred aa a factor to the front of the 
radical. The integrations yield ■ • 
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NATIONAL ADVISORY 
COMMITTEE FOR AERONAUTICS 

Figure 2,- Integration regions for calculating velocity 
potential on surface of thin flat plate v?lng at super- 
sonic speeds. 
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Fleure 6.» Sesims of Integration for calculating velocity 
potential on surface of wing of finite thickness influenced 
by two Independent perturbed flow fields external to wing 
surface . 



